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O , Abstract 

t> ■ 

Q\ ■ The method reducing the solution of the Schrodinger equation for several types of 

power potentials to the solution of the eigenvalue problem for the infinite system of alge- 
Qh! braic equations is developed. The finite truncation of this system provides high accuracy 

results for low-lying levels. The proposed approach is appropriate both for analytic cal- 
culations and for numerical computations. This method allows also to determine the 
spectrum of the Schrodinger-like relativistic equations. The heavy quarkonium (char- 
.^J ' monium and bottomonium) mass spectra for the Cornell potential and the sum of the 

^ ■ Coulomb and oscillator potentials are calculated. The results are in good agreement with 

experimental data. 



1 Introduction 



The solution of the spectral problem for the Schrodinger equation with spherically symmetrical 
potentials is an important problem for the atomic and hadronic spectroscopy. It is well known 
that potential models give a good description of the heavy quarkonium (-0, T) mass spectrum 
0. The interaction potential in such a system is usually assumed to be of the confining type. 
The example is the Cornell potential which contains two terms, one is responsible for the 
Coulomb interaction dominating at small distances and the other corresponds to the string-like 
interaction leading to the confinement. 

The problem of calculating the energy spectrum of the Schrodinger equation with different 
potentials is attracting the attention of physicists for a long time and is being solved by many 
methods: by the direct numerical solution with fixed boundary conditions for wave functions, 
by the variational method and by different modifications combining analytical and numerical 
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approaches. The quasiclassical method works rather well for a wide class of quantum mechanical 
problems |2|, |3[]. However, most of these approaches does not allow to determine explicitly the 
dependence of the solution on the parameters of the potential. 

In this paper we develop further the method proposed in Refs. [|], |5| which accounts for 
the asymptotic behaviour of the wave functions at large distances (or small momenta). This 
approach, in fact, generalizes the method of integral transformations for the kernels of special 
type, determining the correct asymptotic behaviour of the wave functions. As it was shown 
in Refs. 0, [f| such a method is rather efficient for finding the low-lying spectral levels. These 
low-lying levels are of a special interest for the meson and baryon spectroscopy. 

Here we give the exact algorithm for the reduction of the radial Schrodinger equation to 
the algebraic eigenvalue problem, from which it is possible to obtain the spectrum of the initial 
equation using successive approximation method. This approach is rather simple and suitable 
both for analytical calculations and numerical computations and allows to get approximate 
analytical expressions for the low-lying levels. These expressions are suitable for the qualitative 
analysis of the dependence of the spectrum on the parameters of the model potential. It is 
necessary to note that the resulting approximation scheme is rather rapidly convergent. This 
result is of the independent interest and shows the reliability of developed approach. 

In this paper we consider spherically symmetrical confining potentials of the following form 

U k ( r ) = a 2 r k , a > 0, k — 1,2, (1) 

r 

which are important for hadron physics. 

The general construction of the fast converging approximation scheme with the account of 
the algebraic symmetry is given in J|]. Here we apply this method for the potentials of the type 
(§• 

The radial part of the Schrodinger equation can be presented as follows: 



R{r) = XR{r). (2) 



We assume that the total wave function of the solution contains a standard angular part as 
well: 

^nim (r) = — ^ Y lm (6, (f) . (3) 
r 

In this case the effective potential in (||) will be of a well-known type: 

Ueff (r) = l ^±±± + Uk (r) . (4) 

If the potential Uf. (r) does not contain terms more singular at the origin than r -1 (as in 
the case of (|1|)), the behaviour of the radial part of the wave function R (r) ~ r l+1 at the origin 
will be determined by the first term in (d). Let us make the substitution R(r) ~ r l+l F(r), 
where F (0) ^ . Then the equation can be rewritten in the form convenient for further 
transformations 

F" + - {I + 1) F' + [A - U (r)] F = 0. (5) 
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2 The spectral problem for potential U\(r) 

In the case of the potential ([TJ) linearly depending on r the constant a 2 > can be removed 
by the simultaneous change of the variable in eq. (|5]) r — > a~ 2 ^ 3 r and introduction of the 
dimentionless combinations of parameters A* = Aa~ 4,/3 , z* = zaT 2 ! 3 . The solutions of eq. (|5]) 
are sought in the form 



F(r) 



EQ>k h 

k=0 J k 



fk 



2 s 2fc/3 r ' 2fc+4 '+ 5 

3 



41+5 



(6) 



Coefficients fk are chosen from the considerations of symmetrization of the Jacobi matrix, 
constructed from the coefficients a^, and receiving of the correct asymptotical behaviour of the 



wave function ad infinitum R (r) ~ exp 



2 3/2 

3' 



coinciding in this case with the asymptotics 



of the Airy function. Moreover, such a choice of coefficients allows to adjust the kernel of gener- 
alized integral Laplace transformation (Mittag-Leffler transformation). The recursion relation 
obtained by substitution of (^) in (||) is the following: 



/ 2\ /2\ 1 / 3 z* 

{k + 2) a k+2 -{k + 2l + ~) a fe _! + ^-J —-^ B 



'2fc + 4/ + 6 1 



'2 
',3 



-1/3 y- 



(2k+JL±B 1\ 
i\ \ 3 3/ 



g 1 a k+i 



(7) 



Introducing generating function <3> (w) = Y<T=o a-kUJ 2k ^ 3 this relation can be transformed to 
the integro-differential equation: 



(1 - w z J $ w$ 



z* /2\ 4 / 3 , A* 



it(I) 



2 x 2/3 



wzF = 0, 



where 



Kw) 



4,1+1 



dtt— (1 -t) 3 $(wt) 



r (w) 



41+3 



dtt—(l-t) :i $(wt). 



Factoring out explicitly a pole singularity of the function at w 



-41--, 

^(w) = {l-w 2 ) 6 H 



w 



and performing the conformal mapping of the w plane in the unit circle x 

4;+5 » — 41+2 

unknown funcion A (x), where H (u>) = x~s~ (1 — x) 3 A [x), we get an equation: 



5—^ for new 

l+w 



x(l - x) A + 



M + 5 4/ -9 
— 7, 1 7, — % 



A 



2 2/3 z * 

3 4 / 3 r (I 



;i - x y 2/3 1 ( x ) + 



22/332/37; 



^ /*(*). (9) 



Here the integrals /, I* are connected to the function A (x) by the transformations: 



I x) 



I* (x) 



1 - //) n~% (l - -J *A(u)du, 
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1 -u)3u~% (l - -J "~K(u)du. 



According to general theory 0] the necessary asymptotic behaviour of the wave function is 
satisfied if x = is a regular point of the function A (x) . This allows to search for such solutions 
in the form of the nonnegative power series: 



A (x) = Y, A kX k - 



(10) 



k=0 



Performing the integrals /, I* with the help of the direct and inverse integral Mellin transfor- 
mation, it is easy to find the recursion relation for the coefficients A k : 



n + -I + - )A r 
3 6/ 



2 1 \ 00 00 

n - -I + - J A„_i = z F nk A k + A ^ G nk A k . 

6 z/ k=0 k=0 



Here 



1 n 

*nk ~ — y 2^ C n 
n - m=0 



r ( n — m + 



r(m + |) 
r(t) 



r (jfe - m + \ 



r 



r ^/c — m + 1 



nk 



nl 



Hi) 



r(k- 



n+l 



r [n + 



r ( k — ri 



are semi-infinite matrices and A = 5 • 2- 5 / 3 3- 2 / 3 cr 4 / 3 Ar (|) /Y (l),z = 2- 4 / 3 3 2 / 3 a~ 2 / 3 z/Y (|) 
are parameters. The further procedure of getting approximate energy eigenvalues A n consists 
in the finite truncation of the matrices F, G at the dimensions (N + 1) x (N + 1), where 



N = 0, 1,2, 



The resulting N + 1 eigenvalues A r 



0, 1, N correspond to the position 



of the first levels of the spectral problem. The higher is the value of N, the lower is the 
dispersion value. The accuracy of the lowest levels is the highest. In Appendix we illustrate 
the convergence of the proposed approximation by the exactly solvable examples, which are the 
limits of the potential ([I]). There it is shown that the proposed approach exactly reproduces 
the solutions of the Schrodinger equation for the oscilator and the Coulomb potentials. Our 
perturbation results for S wave linear potential rapidly converge to the Airy function zeros. 
For the ground state the 1% accuracy is achieved already for AT = 1. 

For the nonzero z values it is necessary to take higher appproximations in N. However, 
the dependence of the solution on the potential parameters a and z is qualitatively reproduced 
even for N = 0,1, 2. For example, for A^ = the ground state in this potential is 



A 



(iV=0) 
0J 



1.742a 4 / 3 (l + Yj -2.415a 2 / 3 z 



(12) 
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For N = 1 we get two series of levels corresponding to the ground and first exited radial 
states 



a oj —1-\-r, -3 + 1 7, Uoj 



where 



and 



23 43 
— H . 

24 36' 



£> 2 /31, 113 13_\ / 7 , 77 5_\ /l, 11 11W49, 49 7 _ 

4 =q -W + Tu-nV{i2 l + ^-r) + {i2 l + ^-MV 

For N = 2 the solutions of cubic equation correspond to three radial quantum numbers 
0,1,2. 

It is necessary to point out that the reasonable accuracy of these analytical expressions can 
be achieved only in the limited region of parameters z < 2 max (N), a < a max (N). The boundary 
values of parameters are higher for higher values of the matrix truncation parameter N. The 
estimates of the boundary parameter values can be obtained during numerical computations. 
The convergence of this method is guaranteed by the compact nature of the initial differential 
operator and its speed is rather high for carrying numerical computations which do not need 
the extensive computer resources. 



3 The spectral problem for potential U2 if 



Z \* __ X 



As in the previous case, for the potential U2 (r) the substitution r — ► 4^, z* 
removes the coefficient of the highest power of r. Carrying out the similar consideration it is 
easy to get an equation corresponding to Q9) in the form: 



x (1 — x) A 



'l_ 3 _ V 
2 + 4 4 



I 



3_ V 
4 4 



x 



2 v ; 



(14) 



I(x) 



vT* 1 - - 



.T \ 3 



A (u) dtt, 2; 



The recursion relation for the coefficients of the decomposition A (x) = Y^T=o Ak% k is 

l_ 3_A\ . / / 1 A" 
2 + 4 4 



7i + - + t — r I — [n 



4 



^4n.l — Z J! FnkAk, 



k=0 



nk 



2n\ 



£(n+j 



The zero order approximation gives: 

Ar 0) = «(2/ + 3) 



k — n + \ 



-z. 



(15) 



(16) 
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Note that although this approximation for Z = gives the exact result, its accuracy with the 
growth of Z is decreased. Thus for nonzero Z it is necessary to take higher approximations. 
The limiting exactly solvable cases Z — > 0, a — » are considered in Appendix. 

4 Heavy quarkonia mass spectrum 

In this section we apply the developed method of the solution of the Schrodinger equation to 
the calculation of the mass spectra of heavy quarkonia - mesons consisting from heavy quark 
and antiquark (T(bb), ip(cc)). The main assumption of the nonrelativistic quark model 
which is widely used for the heavy hadron description, is the application of the nonrelativistic 
approximation. This approximation works better with the increase of the constituent quark 
mass. The calculations in the framework of relativistic quark models || indicate that non- 
relativistic approximation gives a good description of the static properties of heavy quarkonia 
(such as mass spectrum, radius etc.), while for dynamical properties (such as decays) the ac- 
count of relativistic corrections is considerably more important. In the nonrelativistic quark 
model quarkonium is described by the Schrodinger equation with local spherically symmetrical 
potential, which ensures quark confinement: 

- — Atf(r) + [V(r) - E]^(y) = 0, (17) 

where \x = m 1 m 2 /(m 1 + m 2 ) is the reduced mass for the bound system consisting from the 
quark and atiquark with masses rri\ and m 2 . 

One of the most popular is the Cornell potential which accounts for asymptotic freedom 
in QCD at small distances and the linear rise of the potential with the increasing distance 
between quark and antiquark: 

4 n 

V(r) = — ^ + Ar + B, (18) 

3 r 

where ot s is a strong interaction constant. Sometimes the quadratic confinement is used: 

4 n> 

V(r) = - + Ar 2 + B. (19) 

3 r 

The equation ( |l7|) with potentials (|l8|), ([l9[) reduces to the radial equation @ with the potential 
(U) after substitution: 

g 

\ = 2(i(E - B), a 2 = 2fiA, Z = -a„fi. 

According to our method, the solution of (|T7|) is reduced to the solution of the systems of 
algebraic equations ( |TT1) for the Cornell potential (|18D and (|15|) for potential fll9D . Bottomonium 
(T(bb)) and charmonium (\l/(cc)) mass spectra, obtained by the final truncation of the matrices 
in (|TTD and fljjD at the size 9x9 (N = 8) are given in comparison with experimental data 
in Tables [I], |2|. We use the standard notations for the level centers of gravity: (n + 1)L, where 
n is the radial quantum number. 
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Table 1: Bottomonium (T(bb)) mass spectrum calculated for the sum of linear confining po- 
tential with the Coulomb potential and for the sum of oscillator potential with the Coulomb 
potential (for N = 8). The parameters of the Coulomb+linear potential are A = 0.18 GeV 2 , 
B = —0.29 GeV, a s = 0.39, mj = 4.93 GeV. The parameters of the Coulomb+oscillator 
potential are A = 0.174 GeV 3 , B = -0.05 GeV, a s = 0.345, m h = 4.95 GeV. 



State 


Coulomb+linear 


Coulomb+oscillator 


Experiment 


IS 


9.447 


9.447 


9.4604* 


2S 


10.012 


10.007 


10.023* 


3S 


10.353 


10.389 


10.355* 


AS 


10.629 


10.742 


10.580* 


IP 


9.900 


9.898 


9.900 


2P 


10.260 


10.259 


10.260 


3P 


10.544 


10.593 




ID 


10.155 


10.147 




2D 


10.448 


10.486 





* - 3 5i state. 



We can see from Tables [I], ^ that by fitting the parameters of potentials QT8[ ) and ( ]T9| ) 
it is possible to get a good agreement with the low-lying level masses of charmonium and 
bottomonium for both potentials. Nevertheless, for higher radial excitations the discrepancy 
becomes more pronounced. Unfortunately, these excited levels are near the threshold of the 
open charm and bottom production, which makes the comparison with experimental data more 
complicated. 

The important advantage of the developed method is the possibility to explicitly determine 
the dependence of the spectrum on the parameters of the interquark interaction potential. On 
Figs. 1-3 we plot the dependence of the binding energy (in GeV) on the quark mass tuq (in 
GeV) and the Coulomb potential parameter a s for the Cornell potential flIB|). On Figs. 4-6 we 
present the similar dependences for the sum of the Coulomb and oscillator potentials ([19]). 



5 Mass spectrum of the relativistic Schrodinger-like equa- 
tion 



There are several different equations for relativistic two-particle bound states @-[13|. Most of 
these relativistic equations are the nonlocal integro-differential equations. However, some of 
them can be rationalized and reduced to a local Schrodinger-like equation [14], [0|, which in the 
center of mass frame can be written as follows: 



P" 



2(i R (M) 



+ V(M, r) #(r) = E(M)V(t), 



(20) 
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Table 2: Charmonium (i/j(cc)) mass spectrum calculated for the sum of linear confining potential 
with the Coulomb potential and for the sum of oscillator potential with the Coulomb potential 
(for N — 8). The parameters of the Coulomb+linear potential are A = 0.18 GeV 2 , B = —0.29 
GeV, a s = 0.47, m c = 1.56 GeV. The parameters of the Coulomb+oscillator potential are 
A = 0.174 GeV 3 , B = -0.05 GeV, a s = 0.345, m b = 1.55 GeV. 



State 


Coulomb+linear 


Coulomb+oscillator 


Experiment 


IS 


3.068 


3.070 


3.0675 


2S 


3.697 


3.730 


3.663 


3S 


4.144 


4.331 


4.159* 


IP 


3.526 


3.508 


3.525 


2P 


3.993 


4.095 




3P 


4.383 


4.670 




ID 


3.829 


3.841 


3.770** 


2D 


4.234 


4.415 





* - 3 Si state, ** - 3 L>i state. 



where p 2 = —V 2 = —A; E(M) and /zr(M) are some functions of the bound state mass M; 
potential V(M, r) can also depend on M. 

Equation ( ]2"0"D differes from the ordinary Schrodinger equation ( |TTD by the dependence of 
the reduced mass and potential on energy eigenvalues, which are nontrivial functions of the 
bound state mass. In the case of heavy quarkonia the problem can be simplified by taking the 
nonorelativistic limit of potential in eq. (|20"D as the initial approximation. Solving the obtained 
equation one exactly accounts for relativistic kinematics. The relativistic dynamical effects can 
be taken into account by calculating corrections to the potential using perturbation theory. 



The explicit dependence of the spectrum of the Schrodinger equation (17) on parameters 
obtained here also allows the calculation in the case when the reduced mass fiR and energy 
eigenvalue E depend in the complicated way on the bound state mass M. Thus it is possible 
to solve the relativistic equation of the type fl2HD with the Cornell potential (|TSD . Substituting 
obtained eigenvalues in the initial equation it is possible to reconstruct the wave function and 
to take into account dynamical relativistic corrections using perturbation theory. 

The quasipotential equation [TIJ, which is widely used for the description of bound states 



in quantum field theory (for application to meson properties see [16]]), can be rationalized and 



transformed to the local equation |T5| of the type (^). This equation has the form: 
fb 2 (M) p 2 \ f d 3 q 



\ 2n R 2[i R ) J (2tt 

where the relativistic reduced mass 



*m( P ) = / jAlV(p } q; M)* M (q), (21) 



M 4 - (m 2 - m 2 ) 2 

~ up ' [ ] 
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and the square of relative momentum on the mass shell 



h 2 ( [M 2 - K + m 2 ) 2 ] [M 2 - (mi - m 2 ) 2 ] 

lMJ _ 4M 2 ' 1 j 

V"(p, q; M) is a quasipotential operator. Converting to the coordinate representation it is easy 
to see that quasipotential equation is of the type (|2~0"D with 



F( M \ = b ^Ml = M[M 2 - (mi + m 2 ) 2 ][Af 2 - (m t - m 2 ) 2 ] 
1 J 2/x B 2M[M 4 - (m 2 - m 2 ) 2 ] ' 1 J 



The solution of quasipotential equation (|2lD with the Cornell potential can be found sub- 
stituting relativistic reduced mass ( p2| ) and from ( |2~3| ) in equation (|TT|). Note that for fixed / 



the number of determinant zeros, corresponding to the bound state masses, coincides with the 
matrix order. This is a nontrivial result[] showing that our approach gives the correct number of 
states even for the complicated mass dependence of energy and reduced mass. The calculated 
bottomonium and charmonium masses are given in Tables f| together with the results of nu- 
merical solution of the quasipotential equation with the account of spin-independent relativistic 
corrections of order v 2 /c 2 |I7j and in comparison with experimental data ||. From these tables 



we see that our results for level centers of gravity, calculated using quasipotential equation, 
differ from experimental data by no more than 1% for bottomonium and 3% for charmonium. 
The account of spin-independent relativistic corrections of order v 2 /c 2 [17| further improves the 
agreement with experiment. 



Table 3: Bottomonium (T(bb)) mass spectrum calculated with the quasipotential equation 
for the sum of linear confining potential with the Coulomb potential (for N — 8). NR - 
without account of relativistic corrections to the quasipotential. SI - with the account of spin- 
independent relativistic corrections of order v 2 /c 2 to the quasipotential. Potential parameters 
are A = 0.18 GeV 2 , B = -0.30 GeV, a s = 0.26, m b = 4.88 GeV. 



State 


NR 


SI 


17| 


Experiment 


IS 


9.520 


9.447 


9.4604* 


2S 


9.990 


10.018 


10.023* 


3S 


10.308 


10.349 


10.355* 


AS 


10.575 


10.590 


10.580* 


IP 


9.881 


9.6 


94 


9.900 


2P 


10.211 


10.259 


10.260 


ID 


10.102 


10.157 




2D 


10.386 


10.448 






* _ 


3 Si 


state. 





The resulting equation for bound state masses contains fractional powers. 
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Table 4: Charmonium (ifj(cc)) mass spectrum calculated with the quasipotential equation 
for the sum of linear confining potential with the Coulomb potential (for N — 8). NR - 
without account of relativistic corrections to the quasipotential. SI - with the account of spin- 
independent relativistic corrections of order v 2 /c 2 to the quasipotential. Potential parameters 
are A = 0.18 GeV 2 , B = -0.30 GeV, a s = 0.32, m b = 1.55 GeV. 



State 


NR 


SI 


17 


Experiment 


IS 


3.155 


3.065 


3.0675 


2S 


3.718 


3.669 


3.663 


IP 


3.546 


3.517 


3.525 


ID 


3.833 


3.8 


08 


3.770* 




* _ 


3 ^i 


state. 



6 Conclusions 

In this paper the approach reducing the solution of the Schrodinger equation for power po- 
tentials to the eigenvalue problem for infinite system of algebraic equations is proposed. Such 
potentials are widely used for heavy quarkonium mass spectrum calculations. The finite trun- 
cation (even for small N) of this system provides a high accuracy values for the low-lying levels. 
The explicit dependence of the determinant of the system on potential parameters allows easy 
determination of the solution spectrum dependence on these parameters. It is necessary to 
note that by substituting found eigenvalues in the initial equation it is possible to reconstruct 
the wave function, which has the number of zeros coinciding with the radial excitation number 
(this is in agreement with the known theorem on the wave function zero number). 

We have investigated the particular cases of the studied potentials for which the Schrodinger 
equation has exact solutions. It is shown that for the Coulomb and oscillator potentials the 
developed method reproduces the exact spectrum. For the linear rising potential, in case I = 
(S states), our solutions rapidly converge to the Airy function zeros. 

The mass spectra of charmonium and bottomonium for the Cornell potential (|i~8f) and the 
sum of the Coulomb and oscillator potentials (|T9| ) are calculated on the basis of the developed 
approach. The dependence of spectra on the quark mass uiq and strong interaction coupling 
constant a s is investigated. It is shown that the proposed method allows to determine mass 
spectra of relativistic equations which can be presented in Schrodinger-like form. The mass 
spectra of charmonium and bottomonium for the Cornell potential are calculated on the basis 
of the relativistic quasipotential equation. The obtained results are in accord with experimental 
data. 
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A The limiting cases of the potential ([J) for which the 
Schrodinger equation has exact solutions 

In this appendix we use our method to calculate the spectrum of the Schrodinger equation for 
the limiting cases Z — > or a —>■ of the potential (|l|) and compare the obtained results with 
the exact ones. 

A.l The spectral problem for oscillator potentail 

The oscillator potential is the limiting case Z = of E^O") potential. Setting z = in (|I5D we 
can find the exact solution of this equation, which reproduce the spectrum of the spherically 
symmetrical oscillator 

A n = {An + 21 + 3) a = 2a (V m + ^ , 
where N m is the main quantum number. 



A.2 The spectral problem for linear potential 



The linear potential is the limiting case Z — of Ux(r) potential. Setting also I = (S 
we get the eigenvalues for the spherical linear potential U e ff = r. As it is well 

In Table || we 

ITU) with Z 



states) in (fTl 



known, these solutions should coincide with the zeros of the Airy function [|18 
present the numerical results for the lowest eigenvalues of eq. 



0, / = versus 

the truncation index N. We see that for N values greater than the radial quantum number n 
by 3 the calculated eigenvalues agree with the exact ones within 1% accuracy. This allows to 
conclude that the convergence rate of proposed procedure is rather high. 



A. 3 The spectral problem for the Coulomb potential 

The Coulomb potential is a limiting case a = of potential (0). It is necessary to note that 
the Coulomb potential in Schrodinger equation with potential ([!]) does not give the leading 
asymptotics for r — > oo, which is determined by the confining potential. As a result, we cannot 
get the solution for the Coulomb problem just setting a = in eq. ( [TT]) or eq. (Jl5|). Therefore, 
an additional analysis is necessary. The effective potential in the equation for the radial part 
of the wave function (|2]) consists in this case of two terms of different sign: 



The bound states in the potential (|A.1|) have a negative eigenvalues. Taking this into account, 



we make in eq. (|5|) the substitution r — > r/y |A|, which makes this equation dimensionless, and 
introduce new spectral parameter A = Z/J\X\. After these transformations we get: 



F" + -(l + l)F'+ [-1 + -JF = 0. (A.2) 



r 
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Table 5: The comparison of the S state solutions of eq. (11) for the potential V{r) = r with the 
exact ones determined by the zeros of the Airy {Ai) function, n is a radial quantum number, 
{N + 1) is the truncation index of matrices F and G in eq. (11). 





n 


N 







1 


2 


3 


4 


5 





2 


17747 












1 


2 


33762 


3.90446 










2 


2 


32928 


4.19704 


5.49815 








3 


2 


33298 


4.06720 


5.87930 


7.12835 






4 


2 


33464 


4.08297 


5.46807 


7.57050 


8.84809 




5 


2 


33561 


4.08355 


5.52374 


6.70034 


9.33590 


10.66484 


6 


2 


33623 


4.08466 


5.51523 


6.81654 


7.85466 


11.18931 


7 


2 


33664 


4.08538 


5.51718 


6.77696 


8.03348 


8.98943 


8 


2 


33693 


4.08588 


5.51784 


6.78365 


7.92349 


9.21939 


9 


2 


33714 


4.08626 


5.51828 


6.78389 


7.94408 


8.98263 


10 


2 


33730 


4.08652 


5.51905 


6.78145 


7.95138 


9.01194 


Zeros 
















of Ai 


2.33810 


4.08795 


5.52056 


6.78671 


7.94413 


9.02265 



We seek a solution of eq. ( |A.2 ) in the form of a series in the powers of r: 

A„ 



OO 

n=0 



n + 2l + V 



Substituting this decomposition in eq. ( |A.2j ) we find the recursion relationship: 



(n + l)4»+i - (n + 21 + l)A n _! = -A A n . (A.3) 

We have made the asymptotic symmetrization of the Jacobi matrix in eq. ( A.3Q . Introducing 
the generating function as <&{w) = J2 n A n w n , we get the following equation for this function: 



(l - w 2 ) $' + [A - 2{l + l)w 



$ = 0. 



(A.4) 



The exact solution of eq. ( |A.4| ) is 

= (1 — w) 



-i-i, 



l + w) 



-l-l 



As it follows from general theory, the behaviour of the generating function in the point 
w = 1 for given A determines the asymptotics of the initial problem wave function. Thus only 
for regular function $(«;) in this point, the wave functions are decreasing with the increase of 
r, and hence, they can be normalized in the semi-infinite region. For this purpose the exponent 
should be equal to the integer number n. This condition determines the eigenvalue spectrum: 
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A = 2(n + I + 1). Substituting the explicit expression for the spectral parameter A and taking 
into account the negative sign of the binding energy, it is easy to obtain the well-known Bohr 
formula for the energy of a particle in the Coulomb field: 



Here N m is the main quantum number. 

In this appendix we have shown that our method of the Schrodinger equation solution 
reproduces the exact results for the Coulomb and oscilattor potentials. For S states in linear 
potential our solutions rapidly converge to the zeros of the Airy function. 
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Figure Captions 

Fig. 1 The binding energy as a function of the quark mass and the strong coupling constant 
a s for the Cornell potential. The parameters of the potential are A = 0.18 GeV 2 , B = —0.29 
GeV. Grey patches correspond to S states, patches with grid correspond to P states. 

Fig. 2 The binding energy as a function of the strong coupling constant a s for fixed quark 
mass (rriQ = 4.9 GeV) for the Cornell potential. S states are ploted by solid lines, P states - 
by crosses, D states - by diamonds. 

Fig. 3 The binding energy as a function of the quark mass for fixed strong coupling constant 
(a s = 0.39) for the Cornell potential. S states are ploted by solid lines, P states - by crosses, 
D states - by diamonds. 

Fig. 4 The binding energy as a function of the quark mass and the strong coupling constant 
a s for the sum of the Coulomb and oscillator potentials. The parameters of the potential are 
A = 0.174 GeV 2 , B = —0.05 GeV. Grey patches correspond to S states, patches with grid 
correspond to P states. 

Fig. 5 The binding energy as a function of the strong coupling constant a s for fixed quark 
mass (mQ = 4.9 GeV) for the sum of the Coulomb and oscillator potentials. S states are ploted 
by solid lines, P states - by crosses, D states - by diamonds. 

Fig. 6 The binding energy as a function of the quark mass for fixed strong coupling constant 
(a s = 0.39) for the sum of the Coulomb and oscillator potentials. S states are ploted by solid 
lines, P states - by crosses, D states - by diamonds. 
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Fig. 3 




Fig. 4 
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